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PREFACE 


Part  of  the  Project  RAND  research  program  consists 
of  basic  supporting  studies  in  mathematics.  This 
Memorandum  indicates  a  method  for  treating  the  problem 
of  minimizing  a  type  of  function'  that  frequently  arises 
in  scheduling  and  organization  theory.  The  technique 
presented  here  transforms  the  problem  into  a  dynamic 
programming  one,  which  can  often  be  readily  resolved 
computationally. 
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SUMMARY 


In  a  recent  paper  [1],  L.  Cooper  discusses  a  number 
of  approaches  to  the  problem  of  minimizing  a  function 
such  as 

N 

f(x,y)  -  £  gj  (x  *  y)  • 

i-1  1 


Problems  of  this  nature  arise  frequently  in  connection 
with  scheduling  and  organization  theory. 

The  purpose  of  this  note  is  to  indicate  how  quasi- 


linearization  may  be  used  to  transform  this  into  a 
dynamic  programming  problem  which  in  some  cases  can 
readily  be  resolved  computationally. 
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AN  APPLICATION  OF  DYNAMIC  PROGRAMMING  TO 
LOCATION-ALLOCATION  PROBLEMS 

1.  INTRODUCTION 

In  a  recent  paper,  Cooper  [I]  discusses  a  number  of 
approaches  to  the  problem  of  minimizing  a  function  such  as 

N 

(1.1)  f(x,y)  -  E  gj (x,y) . 

i-1  1 

Problems  of  this  nature  arise  frequently  in  connection 
with  scheduling  and  organization  theory. 

The  purpose  of  this  note  is  to  indicate  how  quasi- 
linearization  may  be  used  to  transform  this  into  a 
dynamic  programming  problem  which  in  some  cases  can 
readily  be  resolved  computationally. 

2.  QUASILINEARIZATION 

To  illustrate  the  method,  let  us  consider  the 
problem  of  minimizing 

(2.1)  f (x,y)  -  E  wt((x  -  +  (y  -  b^2]1/2,  wt  >  0, 

over  all  real  x  and  y.  This  is  a  particular  version  of 
a  location  problem.  (See  [1]  for  a  detailed  discussion 
of  these  questions.)  Other  references  are  also  to  be 
found  there. 


We  begin  with  the  observation  that 


(2.2) 


Hence 


( (x  -  a 


max 

“i^i"1 


t)2  +  (y  -  bt)2)l/2 
I  (x  -  ai)ui  +  <y  -  bi)viJ. 


(2.3)  min  f(x,y)  -  min  max  [E  w .  [  (x  —  a.  )u, 

x,y  x,y  ut,vi  i  1  1  1 

+  (y  -  bi)viJ  1 

■  min  max  [x  E  w.u,  +  y  E  w.v, 
x,y  ui'vi  i  1  i  1  i 


—  (w^a^u^  +  w^b^v^)]. 

# 

Using  calculus ,  we  note  that  the  minimizing  x  and  y 
are  determined  as  solutions  of  the  equations 


(2.4) 


w.  (x  -  a  ) 

E  - 2 - y  y/y  "  0* 

i  ((x-ai)2+  (y  —  bi)z)1/z 


E 

i 


w^(y  -  bi> 

((x  -  a.)2  +  (y 


0. 


0.  Gross  points  out  that  in  certain  £^^^(2^4)  may 
not  hold,  since  the  minimum  point  (x,y)  may  be  at'~'S"p©int 
of  nondifferentiability  of  a  particular  g^(x,y).  For 
example,  the  points  (a^b^)  may  be  (0,0),  (0,1),  (1,0), 
(—1,0),  (0,-1),  with  (0,0)  the  exceptional  point.  This 
difficulty  can  be  avoided  in  a  number  of  ways.  Let  us 
assume  henceforth  that  it  does  not  occur. 


Since  the  values  of 
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u^  and  v^  that  maximize  in  (2.2) 

are  given  by 

_ _ 12.5)  ut  -  (x  -  aj/i  ...  )1/2, 

*•  m 

■  (y  -  bi)/(  •••  )^2, 

we  see  that  (2.4)  implies  that 

(2.6)  £  w.u.  ■  £  w.v.  ■  0 

i  i  1  t  i  i 

for  the  maximizing  and  v^.  (Recall  also  that 

2  2 

+  v£  =  1.)  Hence,  we  may  add  these  conditions  to 
(2.3),  obtaining 

(2.7)  min  f(x,y)  *  min  max 

x,y  x,y  £wiui*£wivi*0 

m 

*  min  max 

x,y  £w^u^»£w^v^*0 

-  max  [-w.  (a.u.  +b.v.)], 

£wiui«=£wivi«0  111  li 

since  x  and  y  have  now  disappeared.  We  still  have  the 
N  constraints 

(2.8)  u?  +  v?  "*  1,  i  -  1,2,  ...,N. 

3.  DYNAMIC  PROGRAMMING 

•» 

To  obtain  a  numerical  procedure  for  maximizing 
N 

(3.1)  -  £  w^(aiu<  +  btvi) 

i=l 


[£[x£w,u,  +  y£w  v. 

i  i  1  i  i  i 

-  wi(aiui  +  bivi)]] 


[-  wi(aiui  +  bivi) ] 


1  *  J'J 


’I - 


subject  to 

the  constraints 

N 

(3.2) (a) 

E  w.u.  -  0, 
i-1  1  1 

N 

(b) 

E  w.v.  -  0, 
i-1  1  1 

(c) 

ui  +  vi  *  1J  1  j 

we  set 


(3.3)  *  cos  0^,  ■  sin  0^, 

so  that  (3.2c)  is  automatically  satisfied.  We  consider 
the  more  general  problem  of  maximizing  (3.1)  subject  to 
the  constraints 

N 

(3.4)  E  w.  cos  0.  -  cw 

i-i  1  11 

N 

E  w.  sin  0.  *  c«. 
i-1  1  1  1 


Setting 


(3.5) 


fN(cl>c2)  ■ 


N 

-*  E  w.(a.  cos  0.  +  b.  sin  0 
i-1  1  1  1  1 


we  obtain  the  recurrence  relation 


(3.6)  fN^cl,c2^  "  max  wN^aN  cos  eN  +  bN  sin  V 


0 


N 


+  wn  cos  ®n*c2  WN 


'N 


sin 


N  >  2,  with  f1(c1,c2) 


immediately  determined. 


Standard  procedures  can  now  be  used  to  determine  the 

* 

sequence  ( f ^(^^2) }  ;  see  [2].  Ultimately,  we  want 
f^(0,0).  There  is  a  small  amount  of  care  to  be  exercised 
in  determining  the  range  of  c^  and  C2  for  each  N. 
This  can  be  determined  inductively. 

4._  UPPER  AND  LOWER  BOUNDS 

Even  if  an  exact  solution  is  not  desired,  the 
transformation  to  the  maximization  problem  of  (2.7)  will 
be  useful  in  furnishing  lower  bounds. 


A  Lagrange  multiplier  can  be  used  to  reduce  the 
problem  to  a  determination  of  a  sequence  of  functions  of 
one  variable. 
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